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transfers. These exploratory computations are carried out in full QCD with 3 flavours of 
sea quarks, but with only two values of m u = which limits our ability to perform the 
chiral extrapolations. The results should therefore be viewed primarily as a demonstration 
of the feasibility of the method. For the K — > ir form factor we compare the new technique 
to the conventional approach and for the pion form factor we assess our results for very 
small momentum transfer with the help of chiral perturbation theory. 
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1. Introduction 

In this paper we investigate the feasibility of using partially twisted boundary conditions 
in lattice computations to evaluate hadronic form factors at any chosen value of the mo- 
mentum transfer. With conventional periodic boundary conditions on the quark fields on 
a spatial lattice of volume 1? the components of the hadrons' momenta are quantized to 
be integer multiples of 2n/L 1 . This leads to a very poor momentum resolution for phe- 
nomenological studies. For example, on a 24 3 lattice with lattice spacing a ~ 0.1 fm, the 
components of the momentum are quantized in steps of about 0.5 GeV . The available mo- 
mentum transfers are therefore also discrete. Using partially twisted boundary conditions 
we show that it is indeed possible to evaluate the form factors at any required momentum 
transfer, provided of course that the hadronic momenta are sufficiently small for lattice 
artefacts to be negligible. We illustrate the technique by evaluating: 

i) the K — > 7r form factors at zero q 2 , where q is the momentum transfer. These 
are required for the determination of the V us matrix element of the CKM matrix. 
In the standard approach to lattice computations of the K — > tt form factor at 
q 2 = 0, first proposed by Becirevic et al. [1, 2] and subsequently employed in a 
number of simulations [3-6], the form factors are calculated (very precisely) at q 2 = 

throughout this paper, L denotes the spatial extent of the lattice. 
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9max (corresponding to the pion and kaon both at rest) and somewhat less precisely 
at other accessible values of q 2 . The results are then interpolated to q 2 = 0. In 
ref. [7], twisted boundary conditions were used in a quenched study to obtain the 
form factors at values of q 2 not accessible with periodic boundary conditions. In 
this paper we demonstrate that by using partially twisted boundary conditions in a 
dynamical simulation it is possible to evaluate the form factors directly at q 2 = 
with comparable errors to the conventional approach but without the need for any 
interpolation in the momentum transfer, thus removing one source of systematic 
error [6]. 

ii) the electromagnetic form factor of the pion at low momentum transfers (from which, 
for example, the charge radius can be determined). In particular we evaluate the form 
factor at values of q 2 below the minimum value obtainable with periodic boundary 
conditions; this minimum is given by \q 2 \ = 2m n (m n — \/m 2 + (2ir/L) 2 ). In contrast 
to recent lattice studies [8-10] this allows therefore for a direct evaluation of the 
charge radius of the pion. 

These are two phenomenologically important examples, but we stress that the techniques 
can also be applied to a wide variety of hadronic matrix elements. The primary aim of this 
paper is to demonstrate the feasibility of the technique. We have therefore used a restricted 
set of quark masses and hence have only a very limited control of the chiral extrapolation. 
Having demonstrated the effectiveness of the technique we will now undertake a large-scale 
computation of the form factors. 

The plan for the remainder of this paper is as follows. In the next section we intro- 
duce all the necessary definitions, correlation functions and ratios of correlation functions 
and then detail the new approaches to compute the scalar and pion form factor. In 
particular, in section 2.3 we explain why it is possible to use a general set of partially 
twisted boundary conditions for the evaluation of the electromagnetic form factor of the 
pion. Section 3 contains the details of our numerical study together with the results and 
finally in section 4 we present a brief summary and outlook. 

2. Description of the Technique 

The matrix element of the vector current between initial and final states consisting of 
pseudoscalar mesons Pi and Pf, respectively, is in general decomposed into two invariant 
form factors: 

<P/(p/)| W(Pi)> = ti iPf (q 2 )(Pi +Pf)n + fp tPf (q 2 )(m-Pf), , (2.1) 

where q = pi — pj is the momentum transfer. For K — > it semileptonic decays is the 
weak current s^f^u, Pi = K and Pf = n, whereas for the electromagnetic form factor of 
the pion is the electromagnetic current, both Pi and Pf are pions and vector current 
conservation implies that /^(q 2 ) = 0. The form factors fp.p f (q 2 ) and fp_ Pf (q 2 ) contain 
the non-perturbative QCD effects. In addition to the matrix elements considered here, 
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form factors for other phenomenologically interesting semileptonic decays, for example for 
B — > 7r, B — > D and D ^ K decays as well as those for decays into a vector final state, 
are also being computed in lattice simulations. As explained in the introduction, with the 
conventional periodic boundary conditions the form factors can only be evaluated at values 
of q 2 such that the components of momenta of the pseudoscalars, pi and pj, are integer 
multiples of 2ir/L. Our aim is to compute the form factors at arbitrary preselected values 
of q 2 and in the following subsections we explain our technique for achieving this. We start 
with a brief introduction to twisted boundary conditions. 

2.1 Twisted Boundary Conditions 

It is well known that the choice of boundary conditions for particles or fields in quantum 
mechanics or field theory in a finite volume governs the momentum spectrum. This obser- 
vation has been exploited in many applications. More recently it has been appreciated that 
by varying the boundary conditions the momentum resolution for lattice phenomenology 
can be significantly improved [7,11-19]. However, if a new simulation (i.e. a new set of 
gauge configurations) were necessary for every choice of momentum, the use of twisted 
boundary conditions would be prohibitively expensive in computing resources and there- 
fore impracticable. In ref. [13] it was demonstrated that for processes without final state 
interactions, such as the form factors studied in this paper, it is sufficient to apply twisted 
boundary conditions only on the valence quarks, whilst using sea quarks defined with 
periodic boundary conditions (see also ref. [14]). In this way the need for new simula- 
tions is avoided and the method becomes practicable. The introduction of such partially 
twisted boundary conditions changes the finite-volume corrections, but, as demonstrated 
in refs. [13,20], they remain exponentially small in the volume and, as is standard, we 
neglect them. 

In our study we use partially twisted boundary conditions, combining gauge field 
configurations generated with sea quarks with periodic boundary conditions with valence 
quarks with twisted boundary conditions, i.e. the valence quarks satisfy 

i;(x k + L) = e ie ^(x k ), (k = 1,2,3) , (2.2) 

where ip is either a strange quark s or a light quark q. By varying 9 we can tune the 
momenta of the mesons continuously. For the purposes of our study it will be sufficient to 
twist only the valence quark in each meson, with the valence antiquark satisfying periodic 
boundary conditions (the generalization to antiquarks with twisted boundary conditions is 
also straightforward). The dispersion relation for the mesons is then 

E = ^m 2 + (p FT + 9/L) 2 , (2.3) 

where m is the mass of the meson and pft is the meson momentum induced by Fourier 
summation (the components of pft are integer multiples of 2n/L). 

For the matrix element in (2.1) with the initial and the final meson carrying momenta 
Pi and pf, respectively, the momentum transfer between the initial and the final state 
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meson is 



Q 2 = (Pi - Pff = { [Ei(Pi) - E f {p f )f - [(p FTji + Oi/L) - (pftj + Of/L)] 2 | . (2.4) 
2.2 Correlation Functions 

In order to determine the form factors we compute two- and three-point correlation func- 
tions. The two-point function is defined by 

Ci(t,pi) = £ e <ft -*< Oi(t, x) Of (o, 0) ) = M (g-^t + e -^(T-t)) , (2 . 5 ) 

where i = 7r or K, the Oj are local pseudoscalar interpolating operators for the corre- 
sponding mesons O n = q^/^q and Ok = s^/^q and we assume that t and T — t (where T 
is the temporal extent of the lattice) are sufficiently large that the correlation function is 
dominated by the lightest state (i.e. the pion or kaon). The constants are given by 
Zi = (Pi\ Oj(0, 0) I 0) . The three-point functions are defined by 

c PiPf (u, t, t f ,ik,p f ) = O f (t f ,x f ) V 4 (t, x) Oj(t u 6) > 

Xf,X 



Z i Z f 



(p f (p f )\v 4 (o)\PM)) 



x{0(t f -t)e- Ei ^- E f^ - 9(t-t 1 )e- E ^ T+t ^- E ^- t ^ , (2.6) 

where Pi j is a pion or a kaon and V4 is the time component of the vector current with 
flavour quantum numbers to allow the P{ — > Pj transition and where we have defined 
Zf = ( I O/(0, 0)| Pj ). Again we assume that all the time intervals are sufficiently large 
for the lightest hadrons to give the dominant contribution. For the remainder of this 
paper we choose to keep ti and tf fixed (with ij = and tf = T/2) and we will therefore 
only explicitly refer to them where necessary 2 and write Cp iPf (t,pi,pf) with just three 
arguments. 

The correctly normalized vector currents in eq. (2.6) are obtained by multiplying 
the local currents used in the numerical simulations by the normalization constant Zy 
defined by 

Zv = 2 cgjtm • (2 - 7> 

where the index B in C^(t, 0,0) implies that the bare vector current is being used. The 
factor of i in eq. (2.7) corresponds to the two terms on the right hand side of eq. (2.5). 



2 In practice we average over the results obtained with different origins of the quark propagators. In 
these cases all the coordinates given here have to be translated accordingly. 
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In order to extract the matrix element (Pf(pf) | V4 (0) | Pi(pi) ) effectively, it is conve- 
nient to define the three ratios: 



Ri,PiP f (Pi,Pf) 



R2,PiP f (Pi,Pf) 



R3,PiP f (Pi,Pf) 

For sufficiently large t and T/2 — i, so that only the lightest mesons contribute significantly 
to each of the correlation functions, each of the three ratios is independent of t and is 
equal to the matrix element ( Pf(pf) | V4 (0) | Pi(pi) )■ Here we are assuming that t is in the 
forward half of the lattice < t < T/2. The correlation functions for t in the backward half, 
T/2 < t < T are readily related to those in the forward half and hence can be combined 
with them to construct the ratios in (2.8). We discuss the quality of the plateaus and the 
numerical determination of the form factors in section 3. 

2.3 The Pion's Form Factor with Twisted Boundary Conditions 

A sketch of the quark-flow diagram for the transition in eq. (2.1), with the final-state meson 
Pf composed of valence quarks (qiq-s) and the initial-state meson with valence quarks (92^3) 
is as follows: 




<73 



For Kf& decays, specifically for the decay K° — » tt + £vi, each of the three valence has a 
different flavour, q\ = u, qi = s and q^ = d, and the partially twisted theory can be readily 
constructed as discussed in ref. [13]. We can therefore introduce three independent twisting 
angles for the three flavours. For the electromagnetic form factor of the pion however, q\ 
has the same flavour as 52, nevertheless it is still possible to use partially twisted boundary 
conditions to evaluate the form factor, with three different twisting angles for the three 
valence quarks, as we now explain 3 . 

(a) We start by imagining that we evaluate the matrix element (vr(pj) | | vr(j)j)) in 
an infinite volume in full QCD with 3 flavours of sea quarks. We assume isospin 
symmetry, and it will be important to note that in this case G-parity implies that 
only the isovector component of the electromagnetic current couples to pions. When 
we consider partially quenched QCD below, this will imply that the vector current is 
composed of valence quark fields. 

3 In the numerical work described in section 3 we choose to keep the twisting angle of the spectator quark 
(#3) equal to zero. 
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/ / Cp rPf (t,pi,p f ) C Pf P t (t,p f ,pi) 

4y/±H*,f \j Cp .(T/2,A)C Pf (T/2,ff f ) » 



f V C, P i P i (*.P'i.Pi)C , p / P / (t,P/,P/)' 



(2. 



Cp i P f (t,j} i ,p f ) j C Pi (T/2-t,p t ) C Pj (t,p f ) C Pf (T/2,p f ) 
C Pf (T/2,p f ) A/ C Pf (T/2-t,p f ) C P% (t,p t ) C Pi (T/2,p,) 



(b) We next consider the partially quenched 3-flavour theory in which m„ = mX = 

= rrij = m]f , where the superscripts V and S label valence and sea, respectively, 
and mf is equal to the physical strange quark mass. The partial quenching arises 
because the mass of the valence strange quark is different to that of the sea strange 
quark. However, since the valence strange quark plays no role in the evaluation of 
(iripf) | Vp | Tr(pi)), this matrix element is correctly given in this theory. 

(c) We now exploit the SU(3) flavour symmetry in the valence sector which implies that, 
for example, 

(ir + (p f ) | *f U | 7T + ( Pi )) = -(TT + ( Pf ) | d-fd\ n + ( Pi )} = (7T + ( Pf ) I u-fs I K°( Pl )) , 

(2.9) 

where the fields in the currents correspond to valence quarks. The pion's form factor 
in this partially quenched theory is therefore equal to that of the K — > it transition 
(the degeneracy of the three flavours of valence quark implies that there is now also 
only a single form factor also for K — > n transitions). 

Up to this point we have not made any approximations. For example, the SU(3) 
flavour symmetry in the valence sector implies that exactly the same diagrams arise in 
chiral perturbation theory (x?T) in the evaluation of the matrix element 
(ir + (pf) | u^s | K°(pi)} in the partially quenched theory as for {ir(pf) \ \ ir(pi)). 
Of course the labelling of the quark content of the mesons may be different; in one 
case we may have a meson with a strange valence quark and in the other an u-quark, 
but as they are degenerate the diagrams give identical contributions (as they must 
by symmetry). In particular, in contrast to the general case for partially quenched 
theories, there are no hairpin contributions here. 

(d) Finally we consider performing the simulations in finite volume, taking all the sea 
quarks to have periodic boundary conditions. For the three flavours of valence quarks 
however, we introduce different twists, 9 U , 9^ and 9 S , which changes the momentum 
spectrum but not the mass spectrum. For propagators in x?T with quark content of 
the form qq (q = u v , d v , s v ), the effects of the twisted boundary condition cancel and 
the spectrum is the same as for the mesons composed of the corresponding sea quarks 
(i.e. as if the boundary conditions were periodic). Thus, in spite of the different 
masses of the valence and sea strange quarks, no double-pole hairpin contributions 
are introduced. 

For valence flavour non-singlet mesons, the momentum spectrum is changed by us- 
ing twisted boundary conditions rather than periodic ones. However, as explicitly 
demonstrated in Appendix A of ref. [13], at one-loop order in the resulting 

summations in finite- volume are equal to the corresponding infinite volume integrals, 
up to exponentially small terms in the volume. Such exponentially small terms are 
also present with periodic boundary conditions and are generally neglected. 

We have therefore established that we can evaluate the pion form factor with 3 different 
twisting angles for the valence quarks (up to the usual exponential precision in the volume). 
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2.4 /^(O) and U„(q 2 ) at Small Momentum Transfers 

We wish to compute the scalar form factor for decays at zero momentum transfer, 
f Kw {0)- The scalar form factor is defined in terms of f^ n and f^- n by: 

fUi 2 ) = fLii 2 ) + ^^JkM\ (2-io) 

and f K M = /+J0). 

Our approach is to use twisted boundary conditions to induce momenta for the pion 
and kaon such that q 2 = 0. A simple way to do this is to take the pion (kaon) to be at rest 
and to tune the momentum of the kaon (pion). We therefore compute the ratios 



R a ,KAPK,0) with \9 K \ = LJ(^^y - m\ and 9 n = 

, 2 \ (2.11) 

and R a , K A0,Pn) with |^| = \^^) 2 - ml and 9 K = 0, 

where a = 1,2,3. The momenta of the mesons are given by px = 9k /L and p n = 9 n /L 
and it can be readily verified that the choices of twisting angles in the two lines of eq.(2.11) 
both correspond to q 2 = 0. 

The required form factor, /^(O), can be obtained directly from a linear combination 
of the ratios in eq. (2.11): 

JkAW- {EK + m ^ ){mK _ K) _ {mK+En){EK _ m ^ {a -1,2,6). (2.12) 

Here Ek (E n ) is the energy of the kaon (pion) corresponding to the momentum induced 
by the twisting angle in the first (second) line of eq.(2.11). 

By using the spatial component of the vector current 14 (k = 1, 2 or 3) it is possible to 
tune the momenta such that q 2 = and qk = so that one obtains the form-factor / + (0) 
directly. We find however, that this procedure leads to a significantly larger statistical 
error. 

The case of the pion's electromagnetic form factor is simpler since current conservation 
implies that f~ w (q 2 ) = so that (dropping the redundant superscript +) 

(vr + ( P/ )|^(0)|7r + ( K )) = f^{q 2 ) ( Pi + Pf ^. (2.13) 

f-Tnriq 2 ) can therefore be directly computed from the ratios in eq.(2.8) by inducing the 
required momenta for the initial- and final-state pions. 

2.5 Comparison with the Conventional Approaches 

One aim of this paper is to compare the precision with which we can determine the form 
factors using the techniques introduced in the preceding subsection with that obtained 
using standard methods. The numerical comparison will be given in section 3, here we 
describe what the conventional approaches are. 
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Recent calculations of /^(O) follow the procedure introduced by Becirevic et al. [1]. 
The scalar form factor at q 2 = is determined using a phenomenologically motivated 
interpolation of /^- 7r (^ 2 ) between the point at q 2 nax = (mx — m^) 2 and points at negative 
values of q 2 which are accessible by Fourier summation with periodic boundary conditions. 
The value of fx^Qmax) 1S readily obtained with excellent statistical precision from 

R 2 -,kA0, 0) = f K Ml„)(m K + m*) . (2.14) 

The interpolation in q 2 (to q 2 = 0) is constrained by computing f^- n (q 2 ) and f^ n (q 2 ) (and 
hence /^(q 2 )) for a variety of values of q 2 . In order to improve precision, where possible 
ratios of correlation functions are used. /^ 7r (? 2 ) and f^^iq 2 ) can be written in the form 

f + K M 2 ) = F( PK , Pn )f K7T (qL x ) ( 1 + l K fl\~l^W )Y , (2-15) 



Er(pk) + E 7T (p w ) 



and 



/^ 2 )=/^(<?W)- (2-16) 

In order to define the quantities F(p,p') and £,(q 2 ) used in these expressions we start by 
defining the ratio of correlation functions 

ft,-. — > \ C k ,Kn(t,PK,Pn) C K K(t,PK,Pn) , . 

RkKPKiPn) = ~ — j—^ — ttt; 77^ — ' 2.17 

Cx7r(i, PK, Pit) C k , KK{t, PK, Pit) 

where k = 1,2,3 is a spatial index and Ck^p i p f {t-,Pi-, P f) is the three point correlation 
function defined in eq. (2.6) with V4 replaced by "54 4 . For time intervals such that only the 
lightest states contribute significantly, R k is independent of t. £ and F are then defined 
by: 



e(? 2 



and 



(E K (Pk) + E K (p n )) (p K +p n ) k + (Ek(Pk) + F 7T (p 7T )) {jp K + P T T ) k R k (PK,Pn) 
(E K (p K ) + E K (p n )) ( PK -Pn) k ~ (E K (p K ) - E n (p n )) (p K +p n ) k Rk(PK,Pwy 

(2.18) 



C Kn (t,p K ,p 7T )C K (t,6)C\(T/2-t,d) Ek(pk) + E n (p n ) 

-^{PK,Pn)- (2.19) 



C Kn (t, 0,0) C K (t,p K ) CV(T/2 - t,p\) m K 



m. 



Having determined /^(q 2 ) and f^^iq 2 ) for a variety of values of q 2 < and at q 2 nax , we 
fit the results to some ansatz for the q 2 behaviour and determine the scalar form factor at 
q 2 = 0. It should be noted that in addition to the systematic uncertainties introduced by 
the choice of ansatz, one is limited to the number of values of q 2 one can use while keeping 
the lattice artefacts small enough for the required precision. 

Recent examples for the use of this approach have been presented in refs. [3-6]. 



4 Since this is the only place where we require the spatial component of the vector current, for simplicity 
of notation we do not introduce the Lorentz index in eq. (2.6). Note however that the minus sign in eq. (2.6) 
becomes a plus sign for spatial indices of the vector current. Cp i p J implicitly corresponds to the matrix 
element of V4. 
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q;6 = 

Figure 1: The three point function Cp i p f (ti,t,tf,Pi,pf) defined in (2.6) in terms of the the quark 
propagator S(t, x\ U, 0) (dashed black) and the generalized quark propagator S'(U, 0; t, x) (solid 
blue). 

3. Results for the Form Factors 

In this section we describe the details of our numerical simulation and present our results. 
3.1 Lattice Parameters 

For the numerical studies presented in this paper we use two ensembles out of the set 
of Nf = 2 + 1 flavour Domain Wall Fermion [21-23] configurations with (L/a) 3 x T/a x 
L s = 16 3 x 32 x 8 which were jointly generated by the UKQCD/RBC collaborations using 
the QCDOC computer [24-27]. A detailed study of the light-hadron spectrum and other 
hadronic quantities using these configurations has recently been reported in ref. [28]. In 
particular, we use the gauge configurations generated with the DBW2 gauge action [29,30] 
at (3 = 0.72. The bare strange quark mass is am s = 0.04 and we use two different ensembles 
with light quark masses ami = 0.02 and ami = 0.01 respectively. The corresponding pion 
and kaon masses are summarized in table 1 and for the inverse lattice spacing we take 
a -1 = 1.6(1) GeV. We use the jackknife technique to estimate the statistical errors. 

We generate the three point functions of type (2.6) by contracting propagators 
S(t, x; ti, 0) from the origin to any point (t, x) with the generalized quark propagator [31] 
defined by 

S'(U, 6;t f ,p f ;t, x) = Y,l5 (S(t, x;t f ,x f h 5 S(t f , x f ;t h 0) e"^) * 75 , (3.1) 

where we suppress the label indicating the twisting angle 5 on the propagators (cf. fig. 1). 
For the present study we generate the generalized propagators in (3.1) with pf = 0, al- 
though it would be straightforward to extend the study to include other values. 

For the form factors in all cases one meson is at rest and the other has a momentum 
which is induced entirely by the twisted boundary conditions. As a result we can evaluate 
all three ratios R a ,Kn(PK,Pn) at a comparable computational cost. A selection of the 
plateaus for the K13 decay is presented for illustration in fig. 2. We find that, for the choice 
of parameters used in this study, ratio Ri,Kn(pK,Pir) and R2,Ktt(pk,Ptt) have the most 

Note that the dagger in the extended propagator reverses the sign of the twisting angle. 



- 9 - 



pronounced plateaus and the fits lead to comparable final statistical errors. The results we 
present in the following have been obtained from fits to Ri,k-k{pk-,P-k)- 

For the electromagnetic form factor of the pion we allow one of the pions to have a 
non-zero value of Fourier momentum (i.e. its momentum is given by (2irri + 9)/L, where ri 
is a vector of integers and 9 is the vector of twisting angles). In order to evaluate the ratios 
Ri,TTinf(PTT t iPwf) arid R2,n i TT f (Pn i ,Pn f ) we would require the generalized propagators in (3.1) 
at non-zero p and here we restrict our computations to the evaluation of R3, 1H -K f {p-K i -,P-w f )- 
A typical result for this ratio is illustrated in fig. 3. 



Ri 



j I - W— 



R 2 



1 I I 1 



Rs 



I I I 



I Zl 1 1 l 1 h- 



*i I I I 1 





T 1 | 


& , 





t/a t/a t/a 

Figure 2: Examples of results for the three different ratios defined in (2.8) for the decay. The 
two rows correspond to the two quark masses am; = 0.01 (upper row) and am; = 0.02 (lower row), 
respectively. In each case the kaon is at rest and the pion has a momentum induced by the twisting 
angle. 




i i i i 



R 3 



t • ■ I 



I I 



t/a t/a 

Figure 3: Example of a typical result for R3 (cf. eq. (2.8)) for the pion form factor for ami = 0.01 
(left) and ami = 0.02 (right). Here the source momentum is \ap\ = 0.39 and the sink momentum is 
0.1. 



3.2 The K -> vr Form Factor, /^(O) 

We computed the i^-form factor using both the conventional and our new approach for 
two values of the light quark mass on 200 configurations with a separation of 10 trajecto- 
ries. The conventional approach has been described in section 2.5 above; for more details 
on the analysis see refs. [1,4,5]. For this conventional calculation, we average all corre- 
lation functions over results from two positions of the propagator source, (0, 0, 0, 0) and 
(8,8,8,16). For the new approach we generate correlation functions using only a single 
source, (0,0,0,0), but averaging over three equivalent twisting angles, i.e. 9i = (#j,0,0), 
(0,#j,0) or (0,0, 9i), with i = it or K. The numerical values of 9 n and 9k obtained using 
eq. (2.11) are presented in table 1. 

In order to choose the twisting angles which correspond to q 2 = using equation 
(2.11) we need to know what the masses of the mesons are. We initially estimated these 
using a subset of 100 configurations from the same ensembles and using only one position 
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ami = 0.02 


ami = 0.01 






0.3765(21) 


0.3002(23) 




arriK 


0.4312(20) 


0.3942(20) 




6 W 


0.838 


1.315 




u K 


0.963 


1.714 


fKw 


a 2 


0.00299(5) 


0.00883(20) 






-0.00012(4) 


0.00016(17) 




le K 


-0.00017(5) 


0.00017(23) 


flTTT 


9tt 


0, 1.6,2.3 


0,1.6,2.3 



Table 1: Summary of kinematical parameters for the studies involving twisted boundary conditions. 



for the source of the propagator. The mean values obtained in this way are amj° = 
0.306, om| 01 = 0.399, am° 02 = 0.374 and am% m = 0.430, to be compared with those 
eventually determined on the full ensemble and averaged over various positions of the 
propagator source given in table 1. The small differences in the central values, together 
with discretization effects in the pion and kaon dispersion relation lead to a small deviation 
in q 2 from 0. We summarize this effect in table 1 where the quoted values for q 2 are obtained 
using eq. (2.4) with the corresponding meson energies determined from fits to the respective 
two point correlation function. 

We present our results in the plots in fig. 4. The left-hand plot shows the data points 
which one obtains from correlation functions with one meson at rest and the other with 
momentum of magnitude \p\ = 0,2ir/L or y/2(2ir/L). Since in addition to f% n (q 2 ) we 
also compute f®x{q 2 ), from equations (2.15)-(2.19) (modified in the obvious way) we obtain 
results for the form factor at four additional values of q 2 for each value of the quark mass. 

The corresponding values of q 2 are determined by using the pion and kaon masses 
obtained from fits to two-point functions as input to the continuum dispersion relation. 
The curves represent a fit using the pole-dominance ansatz 

^ 2 » = T^W- < 3 - 2) 

where M is a parameter fitted from the data. 

The right-hand plot shows a zoom into the region around q 2 = 0. The two data points 
at q 2 > correspond to the results for q 2 iax for which the pion and kaon are both at 
rest; they can be identified by their strikingly small errors. We also display the results for 
fx^iq 2 = 0) obtained from the pole fit (3.2) at each quark mass. In addition the right plot 
of fig. 4 contains the results from the new approach. 

We obtain the following results using the conventional and the new approach: 
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Figure 4: Results for the form factor (ami = 0.02 as full black circles and ami — 0.01 as blue 
squares). Left: All data points entering the conventional approach. Right: Zoom which shows the 
data points for both the new (diamonds) and the conventional approach at q 2 = and the data 
points at <^ ax . 



ami 0.02 0.01 

conventional 0.9951(24) 0.9827(29) 

* K *^ ' new 0.9926(34) 0.9884(34) 

The results for 7x^(0) as determined from the conventional and the new approaches do not 
agree exactly but the discrepancy is statistically not significant. The size of the statistical 
errors is similar in the two approaches, which is an important condition for establishing 
the new technique. The main motivation for the direct approach advocated in this paper 
is to avoid the need for an ansatz with which to perform the q 2 interpolation and this 
is apparently achieved without significantly inflating the error. In addition, it should be 
stressed that the points at negative q 2 are obtained with \p \ = 2ir/L and \p\ = \^22tt/L, 
so that with L/a = 16 one may have concerns about the size of the lattice artefacts at 
these momenta. Our results show that both concerns can now be eliminated by using the 
new approach. 

As stated in the introduction, this is an exploratory study in which we investigate the 
feasibility of the method rather than aim for the ultimate physical results. In particular 
it will be important to check the precision of the direct approach as the mass of the light 
quark is reduced and/or the volume is increased. Reducing the light quark mass leads to an 
increase of Q^ax an d therefore the value for the form factor fK-n-(0) will be more susceptible 
to the choice for the interpolation in q 2 in the conventional approach. Simulating in larger 
volumes, in addition to reducing the finite- volume corrections, enables smaller Fourier 
momenta (with components which are integer multiples of 2ir/L) to be reached, so that 
the data points at non-vanishing momentum in the conventional approach move closer to 
q 2 = and thus better constrain the interpolation. 

3.3 The Electromagnetic Form Factor of the Pion, f^(q 2 ) 

We compute the pion's electromagnetic form factor for two values of the light quark mass 
ami = 0.02 and ami = 0.01 on 200 configurations separated by 10 trajectories in Monte 
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Carlo time. We generate pion three-point functions C nn with all possible mutual com- 
binations of the twisting angles 6\ t 2 as given in table 1 (cf. also fig. 1), with the twist 
being applied only in one direction (i.e. 6^2 = (#1,2)0,0)). After projecting on the Fourier 
momenta of magnitudes |pft| = 0, 2ir/L, and y/2(2iv/L) for the initial state, we are thus 
able to generate data points for the form factor in the entire range from q 2 = to approxi- 
mately lGeV 2 . We average the results obtained for degenerate values of q 2 (see eq. (2.4) ). 
A simplification compared to the scalar form factor is that we do not need Zy] current 
conservation implies that f nn (0) = 1 and this provides the required normalization. 

As was noted in ref. [10], for large values of the initial and/or final pion's momentum 
the argument of the square root in R3 may become negative due to statistical fluctuations 6 . 
To avoid this problem we follow ref. [10] and consider the two-point functions at smaller 
times, making the replacement, 



C n (T/2 - t shift ,p) 
cosh(E 7T (p)t shi{t ) 



(3.3) 



and we find that £ s hift/ a = 7 gives the best results. This removes the occurrence of negative 
arguments in the square root in (2.8) for most of the values of p used in this paper. The 
remaining kinematic points for which we still find negative arguments for the square root 
in the range of time-slices we wish to fit to are removed from our analysis. 

The results are shown in fig. 5. The vertical lines correspond to those values of 
q 2 = 2m- K (m 7T — \/m 2 + (27r/L) 2 ) below which form factors cannot be computed with 
the conventional periodic boundary conditions. The first observation is that our method 



1= 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

-q 2 /GcV 2 



0.1 0.15 

-q 2 /GeV 2 



Figure 5: Left: All the numerical results for the pion form factor (toj = 0.02 circles and m; = 
0.01 squares). Right: Expanded view of low — g 2 -region, now also including experimental results 
(diamonds) [32,33] and / ffff (0) = 1 (star). 



works very well and that we can indeed compute the pion form factor in the low |g| 2 -regime. 

We find that the size of the errors correlates very well with the magnitudes of the 
momenta of the initial and final state pions. This agrees with the observations of ref. [16] 
where the statistical noise as a function of the induced meson momentum was investigated 



6 Since the induced pion and kaon momenta in the determination of the form factor with twisted 
boundary conditions are small this problem does not occur there. 
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numerically in two-point functions and found to increase with the momentum. Similar 
values of q 2 may be obtained by pions with very different momenta, and we find that the 
errors in the form factors are generally smaller if the pions' momenta are smaller. As an 
example consider the two results for the pion form factor at the two neighbouring values 
of -q 2 « 0.14 GeV 2 and 0.16 GeV 2 for ami = 0.02 which we indicate by the dashed circle 
in the right hand plot in fig. 5. In both cases the pion at the sink has a small momentum, 
\ap\ = and 0.1, respectively. At the source however the pion's momentum is 0.25 in the 
one case (smaller error bar) and 0.39 in the other case (larger error bar). In addition, in 
some cases we achieve a reduction in the error by averaging results for the form factor at 
degenerate values of the momentum transfer. 

3.4 Fits to the Pion Form Factor 

Fig. 5 contains the main results of our study for the pion form factor. They are very 
encouraging, clearly demonstrating the feasibility of the method. Ultimately of course, after 
performing our large scale simulation, we will wish to compare our results with experimental 
measurements, but at this stage we can only perform some rudimentary analyses. In 
particular we investigate various fit-ansatze for our present data with the aim of extracting 
the pion's charge radius, 

, df n 
dq 2 

which has also been measured in various experiments. For a qualitative comparison of our 
data to experiment at low values of \q 2 \ we have added the experimental data of refs. [32] 
and [33] to the r.h.s. plot in fig. 5. We note, that there also exist new measurements of the 
Pion form factor at larger values of \q 2 \ [34,35]. 

One approach is to use a pole-dominance (PD) ansatz of the form: 

^ pd(92) = i-« d K)' (3 - 5) 

where the pole mass is related to the pion charge radius by Mp D = 6/(r 2 )pD. Note that 
for lattice data n = 1, whereas for experimental data one either also sets n = 1 or leaves it 
as a free parameter due to uncertainties in the overall normalization (see [36] for example) . 
In order to extract the physical value of the charge radius from lattice data one determines 
Mp D (m,r) for various values of m n and extrapolates the pole mass to the chiral limit (cf. 
e.g. [10]). Since we only have data for two values of the pion mass we will not extrapolate 
Mp D (m T ) here and merely compare the fit results at fixed pion mass with the ones from 
the polynomial ansatze, 

1 1 

fnnM^ 2 ) = 1 + g ( r lhn and /^quad if?) = 1 + g (^)quad <? 2 + C n q A . (3.6) 

Other ansatze are guided by the prediction of the chiral effective theory [37-40] where 
the pion form factor is a well-studied observable [36,41,42]. Here we quote the result at 
next-to-leading order (NLO) [41,42], 

Un,NLo(q 2 ) = 1 + 4 {2L r 9 q 2 + 2H(ml, q 2 ) + H(m 2 K , q 2 )} (3.7) 
Jit 
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where f n = 92.4(5) (3)MeV is the physical pion decay constant [43], Lg the only low energy 
constant relevant at this order of the effective theory and we define 7 

^ (mV) = ^ (s7^ " IW^^^O - 24^ 
f 1 m 2 1 

with i>(m 2 ,q 2 ) = y/q 4 — &q 2 m 2 . 

We have carried out the following fits: 

A) (r 2 ) from fits of (3.5), (3.6) and (3.7) to the lattice data at the unphysical values of 
the quark mass ami = 0.02 and 0.01; 

B) (r 2 ) x and Lg, where the subscript x indicates that the chiral limit has been taken, 
from global fits to the results at the two pion masses using the NLO expression (3.7), 

C) {r 2 ) x from global fits to the lattice results at the two pion masses using the phe- 
nomenologically motivated ansatz 

jvrTr.phcno^ 2 ) = fivn,NLo{q 2 ) + c x q 2 m\ + c 2 q 2 m 2 K + c 3 q 4 . (3.9) 

The results are summarized in fig. 6 and table 2. The plots show how the result of 
each fit changes under a variation of the range in q 2 for the various fit functions discussed 
above. We keep the lowest point in the fit to the q 2 behaviour fixed at q 2 = and vary the 
upper limit. When we quote the lattice results extrapolated to the chiral limit we mean 
the point defined by m w = 140 MeV and uik = 494 MeV. For the final results which we 
quote in table 2 we chose the upper limit of the fit range to be ~ 0.1 GeV 2 in the case of 
the fits A) and B) and about 0.35 GeV 2 for fit C). 

All the fits in case A) for the same quark mass are compatible within errors for small 
values of \q 2 \. The linear fit ansatz starts deviating from the results of the quadratic and 
PD ansatz between 0.2 and 0.3 GeV 2 . The second order polynomial ansatz and the PD 
ansatz turn out to be stable and rather constant in the error over a large range of q 2 . We 
observe a smaller error for the result of the linear fit for small values of the momentum 
transfer and in this region this is therefore the preferred ansatz. 

For case B) we used eq. (3.7) in a global fit over the results for ami = 0.01 and 
ami = 0.02 and try to obtain a prediction for the charge radius in the chiral limit or 
equivalently for the low energy constant Lg (we set fi = m p = 0.77 GeV). For the results 
which we quote we choose values extracted from the fit including all data points in the 
interval q 2 £ [0,0.1]GeV 2 . The results are plotted in the upper r.h.s. plot of fig. 6. The 
blue circles correspond to our lattice results for the charge radius and we indicate the 
PDG value (r 2 ) = 0.452(ll)fm 2 [43] by the horizontal line with error band. Our results 
disagree very significantly with the PDG value indicating that our quark masses are too 

7 Note that the functions TL and v are a slight modification of the functions H. and v in [36]. 



i>(m 2 , q 2 ) log 



2m 2 



q 2 — v(m 2 , q 2 ^ 1 



2m 2 — q 2 + D(m 2 , q 2 ) 



, (3.8) 
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Figure 6: The plots show results from fits to the pion form factor with various fit functions. Left: 

Results for the charge radius at ami = 0.02 (black circles) and ami = 0.01 (blue squares) for linear, 

quadratic and pole dominance ansatz, respectively. Right: Fits B) and C) compared to the PDG 

value for the charge radius. 





fit 


<^>lo.oi/fm 2 


<aio.o 2 /fm 2 


lin. 


eq. (3.6) 


0.294(13) 


0.234(10) 


quad. 


eq. (3.6) 


0.311(16) 


0.242(13) 


PD 


eq. (3.5) 


0.297(28) 


0.239(17) 



fit 


<^>x/fm 2 


L T 9 (m p )\ x 


B) 


0.351(8) 


0.0050(2) 


C) 


0.37(5) 


0.006(2) 



Table 2: Results of pion form factor fits. 



heavy to be described by NLO chiral perturbation theory. This conclusion is reinforced 
by the fact that the values for Lg as determined independently at ami = 0.01 and ami = 
0.02 are significantly different: L 9 | am!= o.oi = 0.0056(3) and Lg\ ami= o.o2 = 0.0048(2). For 
comparison we also quote the value L r 9 = 0.00593(43) determined by Bijnens [36] using 
experimental data together with next-to-next-to-leading order chiral perturabation theory. 

The results for fit C) are summarized in the lower r.h.s. plot in fig. 6. This global fit of 
the phenomenologically motivated ansatz (3.9) leads to results for the charge radius with 
larger errors. The results which we quote in table 2 correspond to the values extracted 
from the fit including all data points in the interval q 2 G [0, 0.35]GeV 2 . 
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4. Conclusions 



In this paper we have demonstrated the feasibility of using partially twisted boundary 
conditions to compute weak and electromagnetic form factors. The technique allows the 
form factors to be evaluated for any choice of momentum transfer, q, within the range of 
hadronic momenta such that lattice artefacts are small. We have illustrated the techniques 
for two phenomenologically interesting quantities, the evaluation of the form factors for 
Kt 3 decays directly at q 2 = (thus avoiding the need for an interpolation in q 2 ) and of the 
electromagnetic form factor of the pion at small momentum transfers (thus enabling us to 
evaluate the pion's charge radius without extrapolation of the results from large values of 
-q 2 )- 

The computations described here were intended to demonstrate the proof of concept, 
and as such were limited both in statistics and in the lattice parameters. In particular we 
only compute the form factors for two values of m u ^ and hence any investigation of the 
chiral behaviour is restricted. The next step will be to implement the technique in a large 
scale simulation which will enable us to overcome the limitations of this feasibility study 
and to obtain results with the systematic uncertainties under control. 
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